Dynamics of excitations in a one-dimensional Bose liquid 
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We show that the dynamic structure factor of a one-dimensional Bose hquid has a power-law 
singularity defining the main mode of collective excitations. Using the Lieb-Liniger model, we 
evaluate the corresponding exponent as a function of the wave vector and the interaction strength. 
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Progress in the ability to manipulate ultracold atomic 
gases stimulates the interest in fundamental properties 
of one-dimensional (ID) Bose liquids The quan- 

tity characterising the collective excitations in these sys- 
tems, the dynamic structure factor (DSF), is now di- 
rectly accessible experimentally using the Bragg spec- 
troscopy technique [2]. Already the very first such mea- 
surements 3] clearly showed that the resonance in DSF 
is wider in ID than it is in higher dimensions. The goal 
of this Letter is to elucidate the nature of the resonance 
in a ID system of interacting bosons. 

In the absence of interactions, bosons occupy the 
lowest-energy single-particle state at zero temperature. 
An external field that couples to the particle density 
would excite bosons from the ground state. The cor- 
responding absorption spectrum reflects the free boson's 
dispersion relation e(g). Accordingly, DSF at zero tem- 
perature is given by S{q,Lj) oc 6(u) — e(g)). 

In dimensions higher than one, this behavior remains 
largely intact even in the presence of interactions. Bosons 
still form a condensate, and excitations of the system 
are very well described in terms of Bogoliubov quasi- 
particles Interactions merely affect their spectrum: 
e{q) cc q at small q. The quasiparticle decay rate scales 
with q as I/t^ oc [^]; hence, the quasiparticle peak in 
S{q,LLi) at w = e{q) is well defined, l/r^ ^ e(q). 

In ID, the effect of interactions is dramatic: quantum 
fluctuations destroy the condensate. Long-wavelength 
{q 0) excitations of a ID Bose liquid are often de- 
scribed in hydrodynamic approximation Q (see for a 
recent review). However, the shape of the peak in DSF 
can not be addressed using this approach: in hydrody- 
namics the peak has zero width. 

In this Letter we study DSF of a ID Bose liquid be- 
yond the hydrodynamic approximation. We consider the 
Lieb-Liniger (LL) model 7]: N identical spinless bosons 
with contact repulsive interaction placed on a ring with 
circumference L, 
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The model is integrable 0, Q • The integrability allows 



one to relate the parameters of the hydrodynamic de- 
scription Q , the sound velocity v and the parameter 
K, to the concentration n = N/L and the dimensionless 
interaction strength 7 = mc/n [9|. Finding dynamic cor- 
relation functions, such as DSF, in a closed form remains 
a challenge The most impressive progress so far was 
achieved by combining a finite- TV numerics with the Al- 
gebraic Bethe Ansatz J_0,]. Here we study the singular 
behavior of DSF analytically. 
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FIG. 1: (a) Shaded area indicates the region in (Lj,g)-plane 
where S{q, tj) 7^ at zero temperature. DSF exhibits power- 
law singularities along the solid lines, (b) Sketch of the de- 
pendence of the structure factor on u at a fixed q < 2nn. 

DSF is defined by 

S{q,ij)^ y"dxdie'('^*-«^) (p(x,t)p(0,0)), (2) 

where p{x) — J2i ~ ^i) the density operator. We 
show that DSF exhibits power-law singularities at the 
Lieb's modes £1,2(9) 0, EHj see Fig. [TJa). In particular, 
DSF diverges at w ^ ^lil) 



S{q,uj) ~ — 

q 
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ei 



[e{ei-uj) + i^ie{u;-ei)], (3) 



see Fig. [TJb). Here 5e{q) = minjci — e2,vq}. Note that 
the divergence occurs within the continuum. 

The exponent fii and the coefficient vi in Eq. ^ de- 
pend on the dimensionless momentum Q — q/mc and the 
interaction strength 7. We were able to compute jji and 
h'l in two limiting cases: 
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for Q ^ {'^K) ^ niax{l,7 ^/^} and arbitrary 7, and 

sin((5747r) 
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for 7 1 and arbitrary Q (here 5 = 2 arctan Q) . Ac- 
cording to Eq. ^1 « 2(5/7r at Q — > 0; we expect that 
/ii cx Q at small Q for any 7. 

The line oj = £1(9) has a "replica", w = £2(9), at 5 > 
27rri. Here DSF does not diverge, but still has a power- 
law non-analyticity. The singular part of DSF has the 
form dS{q,Lo) oc \uj — t2\^^ with the exponent 



/i2 = 2K+{2K)-^ - 1 
at arbitrary 7 and Q ^ max{l,7^^}, and 
fi2^{S/n){l + S/27r) 
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at 7 ^ 1 and arbitrary Q. 

Lieb's hole-like (according to Bethe-ansatz classifica- 
tion) mode [m £2(9) serves as the lower boundary of the 
support of S{q,uj) at q < 27rn, see Fig. [H Here DSF is 
given by 



S{q, 
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For 7^1 the exponent ^2 here is given by Eq. ([7]). 

Equations jS])-® represent the main result of this Let- 
ter. The shape of S{q, lj) near uj = £i('?), see Eq. ([3]), dif- 
fers qualitatively from the Lorentzian quasiparticle peak 
in higher dimensions. In ID, the collective mode is char- 
acterized by a power-law divergence of S{q,uj). This 
divergence is protected by the integrability and associ- 
ated with it absense of three-particle collisions @. It is 
smeared only at a finite temperature T, 



max 
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An apparent saturation of the height of the peak with the 
decrease of T would provide a direct measure of three- 
particle scattering (absent in LL model) or recombina- 
tion [l3| rates. 

In the remainder of the Letter we outline the derivation 
of the above results. We start with the limit of large q. 
Consider the state |0), where = J^k'^l+q'^k 
Fourier component of the density operator {ipj, creates a 
boson with momentum p), and |0) is the ground state of 
the Bose liquid. Without interactions, all bosons in |0) 
occupy the single-particle state with fc = 0. The operator 
annihilates one such boson while creating another in 
the empty state with momentum q. With interactions 
present, the occupation number falls off [ist rapidly with 
k at k > mv. Therefore, for q ^ mv the state pj |0) still 
contains a single particle at momentum close to q, as well 
as a "hole" in the quasi-condensate with much smaller 
momentum. This observation suggests to approximate 



pj « / dx d\x)ip{x), 



(10) 



where d^{x) = L^^^'^J2\k\<ko^~^'''' '^l+k '^'"^^.tes a high- 
momentum particle and 4^{x) creates a long- wavelength 
hole; here, ko ~ mv is the high-momentum cutoff. The 
d-particle is described by the Hamiltonian 

Hd= dxd'^{x)[ei{q) - ivdda;]d{x), Vd = q/m. (11) 



Here we took into account that £i(p) « p'^/2m at large 



111 ] and linearized the dispersion relation around p = q. 



We treat the long-wavelength bosons in the conventional 
hydrodynamic approximation [^, Q , 



(12) 



The fields (p, t9 obey [(p{x), ^{y)] = i{'^/2) sgn(a; — y) and 
their dynamics is governed by the Hamiltonian 
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Eqs. ©, (Uni), and in]) yield DSF in the form 

S{q,uj)= [dxdte''^'{B{x,t)Bl {0,0)) (14) 



with B^x) oc d'f(a;)e^''(^). Evaluation of Eq. ^ with 
the quadratic Hamiltonian + Hq is straightforward 
and yields Eq. ([3]) with /ii and i^i given by Eq. Q. The 
decomposition Eq. pUj) is applicable for q :$> ko ^ mv, 
hence the restriction on Q in Eq. (|4]). On the other hand, 
the constraint |fc| < kg on the momentum of d-particle 
limits the applicability of Eq. ([3|) to |a; — £i(g)| < qv. 

We now extend the above derivation to the vicinity 
of the mode £2(9) at q ^ max{27rn, 717}. At these mo- 
menta, £2(9) = £1(9— 27rri) is areplicaof mode £i(g). At a 
given energy w w £2(9) the relevant excitation includes, in 
addition to d-particle, the 27rn-momentum excitation of 
the quasi-condensate Eq. (fT3|) . In hydrodynamics [1, 0] , 
such excitation corresponds to ip{x) cx 
instead of Eq. DSF is stiU given by Eq. ^ with 

St(.x) oc dt(a;)e*''(^)-2«[™^+'^(^)l and with the replace- 
ment £1 £2 in Eq. (fTTIl . Evaluation of Eq. (fT4|) then 
yields a power law for the singular part of S{q,uj) with 
the exponent fi2 given by Eq. 

At 7 ^ 1, one can take an advantage of the exact 
mapping [l3| of the LL model onto fermions with 6"{x) 
interaction. The mapping generalizes the famous dual- 
ity [isl between impenetrable bosons and free fermions 
and is based on the elementary identity 



2 arctan (p/mc) = tt — 2 arctan (mc/p). 



(15) 



The second term in the r.h.s. here is the scattering 
phase shift 9 sip) of the symmetric wave function [3] of 
two particles with relative momentum p interacting via 
Vb = cS{x) potential. Adding tt to 6s converts the sym- 
metric wave function into the antisymmetric one. On the 
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other hand, the l.h.s. of ITSl) is the phase shift Oa{p) of 
the antisymmetric wave function of two particles inter- 
acting via potential Vp = —2/{m?c)5"{x). In view of 
the integrability of the LL model, the two-particle phase 
shifts contain a complete information about the Bethe 
ansatz wave function of the many-body problem. Thus, 
for any bosonic eigenstate of the LL model ^ there is a 
dual fermionic eigenstate of the Hamiltonian 



Hp 



9m ' ^ rin^"^ n^'^r 



2m ^ dx - 

j = l J 



5"{x^-x,) (16) 



that has the same energy. The two wave functions co- 
incide in one of the sectors, say xi < X2 ■ ■ ■ < xm, but 
differ by their symmetry with respect to the permuta- 
tion of particles' coordinates. Since the density operator 
does not permute particles, its matrix elements between 
any two many-body eigenstates of Eq. ([T]) are identical to 
those evaluated with the corresponding dual eigenstates 
of Hp. In particular, the DSF for the LL model coincides 
with that for the fermionic model Eq. (fT6|) . 

It is convenient to rewrite Eq. in the second- 

quantized representation, 
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Vk 
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Here the operator creates a fermion with momentum 

p and energy = p^/2m and pk = Y.p V'p-fcV'p- 

Strong repulsion between the original bosons corre- 
sponds to a weak interaction in the dual fermionic model 
Eqs. (dH), ini). In the limit c ^ oo Eqs. HI]), HH) 
describe free fermions. In this limit the structure factor 
differs from zero only in a finite interval, e2 < w < ei with 
£1,2(9) — voq ± q^/2m. A weak (oc I/7) residual inter- 
action between fermions leads to corrections to ei.2; for 
example, the Fermi velocity vq is replaced by the sound 
velocity v ^ vq/K ^ uo(l — 4/7) [l,|iH. Rather than dis- 
cussing these modifications, we concentrate here on the 
singularities in S{q,uj). 

DSF is proportional to the dissipative response to a 
field that couples to density. In the fermionic represen- 
tation, the absorption of a quantum with energy to and 
momentum q is due to excitation of particle-hole pairs; 
there is just one such pair in the limit c— >oo. Atw^ei, 
the hole is created just below the Fermi level while the 
particle has momentum close to kp + q; here kp — Trn/m 
is the Fermi momentum. In the presence of interactions, 
such process is accompanied by a creation of multiple 
low-energy particle-hole pairs near the two Fermi points 
p — ±kp. Similar to the well-known phenomenon of the 
Fermi ed ge s ingularity in the X-ray absorption spectra 
of metals [l7| , the proliferation of low-energy pairs leads 
to power-law singularities in the response function at the 
edges of the spectral support. 



The Fermi edge singularity relies crucially on the 
sharpness of the distribution function which is smeared 
at a finite temperature, hence Eq. ([9]). Note that in an 
integrable model Eq. ([T6|) there is no relaxation of ex- 
cited fermions [ISj as three-particle collisions are absent. 
Therefore, at T = there is no smearing of the power-law 
singularities in S{q, uj) even at finite uj [is'. 

The derivation of Eqs. dS]), dS]), d?]), and dH]) follows the 
method of Ref. [l^l- Consider first the limit w ei. We 
truncate the continuum of single-particle states to three 
narrow subbands of the width ko -C q: d-subband 
around p = q that hosts a single particle in the final state 
of the transition, and two subbands, a = ±, around the 
right/left Fermi points p = ikp that accommodate low- 
energy particle-hole pairs. After linearization of the spec- 
trum within each subband, the resulting effective Hamil- 
tonian takes the form 



H — Hd + Ho + Hint 
with Hd given by Eq. (1111) and 



(18) 



Ho = Jdx E.V'i(^)[- 



with iPaix) = L-i/2^^e*(fe-"feF)2 
the r.h.s. of Eq. dlS]) describes interaction. 

Hint = Ug J dx Pa{x)pd{x), 



-iavda:]lpa{x) (19) 

ij^k- The last term in 



(20) 



where pa = V'iV'a ^^'^ coupling constants Ua are 
related to the parameters of the initial Hamiltonian (|17|). 
see Eq. dSIl) below. 

Note that Hint does not include the direct interaction 
between the right and left movers. Indeed, the corre- 
sponding coupling constant V2kp = 871110/7 is small in 
the limit 7 3> 1. In the absence of such interaction, the 
remaining coupling constants Ua are set by the require- 
ment [1^ that the two-particle scattering phase shifts for 
the effective Hamiltonian (|18p - d20)) with linearized spec- 
trum reproduce those for the original model dlB]) . (|17|) . 
In the latter case, the phase shifts 6± = 9a{q + kp T kp) 
(see [3]) are given by 



S± = 2 arctan 



2{vd T v) 



In the limit 7 3> 1 taken at a constant Q = q/mc^ one 
finds 5± = 5. In order to reproduce these phase shifts, 
the coupling constants in Eq. dlO]) must be equal to 



[vd 6 — 2 arctan Q. 



(21) 



In terms of the effective Hamiltonian (fT8)) - d2T|) . the struc- 
ture factor is given byEq. dH with fit (3;) ^ d'' {x)il:+{x). 
Following the steps familiar from the theory of the Fermi- 
edge singularity [2l|, we arrive at Eq. dSl) with the expo- 
nent given by Eq. 
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The power law Eq. dH) with exponent of Eq. ([7]) is ob- 
tained in a similar fashion. The only difference is that 
at q < 2kp the d-subband is centered at momentum 
p = kp — q (i.e., it is below the Fermi level) and hosts a 
single hole with velocity Vd = v — q/m. The hole is rel- 
atively slow, \vd\ < V, which leads [2^ to S{q,uj) — 
at a; < £2(9)- With q 2kp — 0, the center of d- 
subband is approaching —kp. At larger one returns 
to the particle-like d-subband, but the density operator 
fit in Eq. (mi) IS now given 

The "bosonic" route of evaluation of DSF described 
first, and the "fermionic" one described second, have a 
common region of applicability corresponding to both Q 
and 7 being large. In this limit the fermionic calculation 
yields fii -> 1/2 and /i2 ^ 3/2, see Eqs. jS]), ©. This 
is in agreement with the strong- repulsion limit {K 1) 
of the result of the bosonic calculation, see Eqs. ([4]), ([6]). 
Note that at the special point q = 2Tm the exponent ^2 
can be found using the hydrodynamic approximation (in- 
deed, £2 ^ at (? — > 2?™, hence oj — > £2 limit is accessible 
within the effective low-energy description). The hydro- 
dynamics yields ^2 = K — iKiA/^^ , in agreement 
with the corresponding limit of Eq. ([7]) . 

The peculiarity of ID Bose liquid is that an arbitrar- 
ily weak repulsion between particles destroys condensa- 
tion. This renders the perturbation theory developed for 
higher dimensions [3] inapplicable. The well-known alter- 
native method based on the hydrodynamic description of 
low-energy excitations [S] also has its limitations, yield- 
ing infinitely narrow resonance in DSF at small g (6|. In 
this Letter we demonstrated the existence of power-law 
singularities in DSF. The two complementary methods of 
analytic evaluation of DSF developed here allowed to find 
the corresponding exponents in several regimes. Evalu- 
ation of the exponents in the entire range of parameters 
remains a challenging problem. 
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